Surface magnetic ordering in topological insulators with bulk magnetic dopants 
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We show that a three dimensional topological insulator doped with magnetic impurities in the 
bulk can have a regime where the surface is magnetically ordered but the bulk is not. This is in 
contrast to conventional materials where bulk ordered phases are typically more robust than surface 
ordered phases. The difference originates from the topologically protected gap less surface states 
characteristic of topological insulators. We study the problem using a mean field approach in two 
concrete models that give the same qualitative result, with some interesting differences. Our findings 
could help explain recent experimental results showing the emergence of a spectral gap in the surface 
state of Bi2Se3 doped with Mn or Fe atoms, but with no measurable bulk magnetism. 



I. INTRODUCTION 

In recent years, the field of topological insulators (TI) 
has attracted much attention and research in the con- 
densed matter community.'^l'^ The advance has been 
rapid, on both the theoretical and experimental fronts, 
however, many challenges still remain. Among these 
perhaps the most important is gaining experimental 
control over the bulk and surface conduction in three- 
dimensional TIs. In this paper we address one aspect 
of this challenge that is associated with a magnetically 
induced excitation gap in the topologically protected sur- 
face states. 

Topological insulators are bulk insulators in 2D or 3D 
with strong spin orbit coupling (SOC) and protected 
gapless surface states.l^H^ The topological protection of 
the surface states arises due to time reversal invari- 
ance (TRI). The surface states are conducting, have a 
characteristic linear (Dirac) dispersion and exhibit spin- 
momentum locking. The most studied and most promis- 
ing 3D TI is the semiconducting thermoelectric Bi2Se3 , 
with a relatively large band gap of ~ 0.3eV and a simple 
surface state consisting of a single Dirac cone.^^ The 
spectrum of Bi2Se3and other TIs has been studied us- 
ing angle resolved spectroscomB^^^ (ARPES) and scan- 
ning tunnelling microscopy,^^ showing that the sur- 
face states form an almost ideal Dirac cone, illustrated 
in Fig. [l^, familiar from studies of graphene.^^ 

Breaking TRI, for example by adding magnetic 
dopants (as we shall discuss), is expected to open a gap in 
the spectrum of the surface states. The resulting spec- 
trum then resembles that of a 'massive' Dirac fermion 
(Fig. There is considerable interest in having a sys- 
tem with an odd number of massive Dirac fermions, since 
it is predicted to exhibit many interesting topological 
phenomena, including the half quantum Hall effect on 
the surface {e^ /2h Hall conductance),^ the image mag- 
netic monopole (an electric charge adjacent to a TI re- 
sults in the field of a magnetic monopole embedded in 
the TI),^^^ and a Kerr/Faraday angle quantization in 
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FIG. 1: (a) Energy spectrum of a massless Dirac fermion. The 
bottom (red) cone is the valence band, fully occupied at half 
filling. The top (green) cone is the conduction band, which 
is assumed to be vacant, (b) Massive Dirac fermion. This 
figure shows how the opening of a gap tends to lower the free 
energy. 



units of the fine structure constant A tunable gap 
would also allow the control of the surface transport, 
and could in addition lead to unique practical applica- 
tions associated with purely electric control of the surface 
magnetization.^^ 

A signature of the massive Dirac fermion has been 
observed recently using ARPES in magnetically doped 
Bi2Se3 although the interesting effects associated 

with it have yet to be seen in a laboratory. A surprising 
feature of these experiments is that the gap in the surface 
spectrum appears without bulk magnetic ordering, even 
though the dopants are uniformly distributed everywhere 
in the 3D sample. These findings raise several important 
questions concerning the precise conditions under which 
TRI-breaking perturbations open up a gap. Can a gap 
open in the surface state of a TI in a TRI-broken phase 
which however lacks global magnetic ordering? Although 
we know of no systematic study of this problem, simple 
arguments suggest that unordered magnetic moments do 
not open a gap. Consider creating such a disordered state 
from a uniform 2D ferromagnet (EM) in the surface of a 
TI by introducing domains with opposite magnetization 
(taken to point in the direction perpendicular to the sur- 
face). It is well known that the resulting domain walls 



2 



carry topologically protected gapless fermionic modes.^ 
As the number of the domains grows so does the density 
of the low-energy fermion modes, ultimately presumably 
recovering the 2D gapless state characteristic of the sys- 
tem with unbroken TRI. The above argument thus sug- 
gests that uniform magnetic ordering over large domains 
is necessary to gap out the surface modes in a TL 

In this paper, we lay out the hypothesis that a tem- 
perature window exists in which the surface of a mag- 
netically doped TI is magnetically ordered but the bulk 
is not. We present a simple and intuitive argument why 
this is so, and back it up by a mean- field calculation 
for two simple tight binding TI models: a cubic-lattice 
regularized Bi2Se3 and a model on the perovskite lattice. 
Our results show that indeed a sizeable regime such as 
described above could exist in real TIs, and this indi- 
cates a possibl e phy sical explanation for the results seen 
in experiments. ^^^'^ 



II. SURFACE MAGNETIC ORDERING IN 
TOPOLOGICAL INSULATORS 

A. Surface doping 

The most natural way to attempt to open up a gap 
in the surface state of a TI is to coat the surface with 
a ferromagnetic material, with magnetization perpendic- 
ular to the surface. Theoretically, this causes a gap to 
open u p, pro portional to the magnetization of the FM 
coating .f^SI^ To illustrate this point, consider the effec- 
tive low energy Hamiltonian for electrons on the surface 
of a 3D TI that lies parallel to the x-y plane 



(1) 



where v is the Fermi velocity, and di are Pauli matrices 
in the spin subspace. If we coat this surface with a ferro- 
magnetic coating with magnetization M = Mz then we 
get an additional term in the Hamiltonian 



H = Ho + JM 
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where J is the exchange coupling strength. Since H is a 
sum of anti-commuting matrices we can write the spec- 
trum down immediately 



Ek = ±^v^k^ + (JM/2)2 



(3) 



where k'^ = kl ^ k^. 

We see that a gap of size JM has opened up. However, 
even though from a theoretical standpoint this proposi- 
tion seems promising, experimentally it has proven very 
difficult to fabricate a sample with the requisite prop- 
erties. Two key challenges need to be overcome: first, 
to observe most of the interesting surface phenomena, 
one requires the surface to remain insulating; however 
most ferromagnets in nature are metallic. Second, for a 
ferromagnet in a thin-film geometry, the magnetization 



vector usually lies in the plane, whereas a perpendicular 
magnetization is required to open up a gap in the TI sur- 
face state. To the best of our knowledge, this has yet to 
be achieved in an experiment, although some theoretical 
work has been done on this topic 



B. Bulk Doping 

If surface doping with magnetic impurities fails, i t is 
natural to try bulk doping. In ARPES experiment s^^'^ 
it was found that doping the bulk with non-magnetic im- 
purities (such as Ca, Sn and TI) did not result in a gap 
in the Dirac cone, as expected since they do not break 
TRI. Conversely, doping with magnetic impurities, for 
example Bi2-a;Fea;Se3 , resulted in a spectral gap that in- 
creased with the concentration of magnetic dopants x, 
with a gap of 60meV for x = 0.25 (the bulk gap for 
Bi2Se3 is ~ 0.3eV). For the magnetic dopants Fe and Mn 
it was found that, at least for small x, the bulk was para- 
magnetic, while for the undoped samples the bulk was 
found to be diamagnetic. The magnetization measure- 
ments were not sensitive to the surface. 

This raises the question of magnetic ordering in the 
bulk versus the surface. In general, ordered phenomena 
in lower dimensions are more fragile (T^^ > T^^), for 
example, the XY model and Heisenberg models - in ID 
they do not order at any temperature, in 2D they order 
only for T = and in 3D they order for T < Tc. This 
is also the case for superconducting order and general 
stability of lattices. However, in the case of a TI, we will 
argue that it is possible that T^""^^ < T^'^'^ Therefore, 
there is a regime T^^^^ < T < T^^^^ in which the bulk is 
unordered (paramagnetic) and the surface is ordered (for 
example, ferromagnetic). 

To illustrate why this could be the case, first recall that 
magnetic ordering with the magnetization perpendicular 
to the TI surface implies opening of a gap in the spec- 
trum of the surface states. This can be seen directly 
from Eqs. (2|3). Now consider the ungapped surface 
spectrum, assuming half filling, so that the surface va- 
lence band is fully occupied and the conduction band is 
empty (Fig. [l^). Gapping the surface states causes the 
occupied states to move down in energy (Fig.[l]3), so the 
total kinetic energy decreases. Therefore, the formation 
of a surface gap is favourable. If the chemical potential 
is shifted either up or down then the net gain in kinetic 
energy is diminished and we expect T^^^^ to decrease. 

Contrast this with the situation in the insulating bulk, 
which is gapped to begin with. In an ordinary insulator 
with negligible spin-orbit coupling it is not possible to 
generate magnetization in the initially spin-degenerate 
bands without first closing the gap. Equivalent ly, one 
may recall that the spin susceptibility of an ordinary 
insulator with a negligible spin-orbit coupling vanishes. 
In the bulk of a topological insulator the situation is 
more complicated as a result of the strong spin-orbit cou- 
pling that is necessary for the occurrence of the topo- 
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logical phase. In this case, magnetic susceptibility can 
be significant,^^ and can lead to bulk magnetic states at 
non-zero temperatures. Nevertheless, in this study we 
find that generically the surface critical temperature for 
the formation of magnetic order exceeds the bulk critical 
temperature. 

Equivalently, we can imagine integrating out the elec- 
trons, and consider that the range of the RKKY interac- 
tion between the dopants is inversely proportional to the 
gap. Since the gap is small for the surface and large for 
the bulk, the range of the interaction is much longer for 
the surface, and hence long range order is expected to be 
stronger on the surface. 



III. FORMALISM 

A. Mean field theory for the bulk system 

We consider first just the 3D bulk of the material, 
and study the coupling of electrons to magnetic dopants 
within a mean field approximation. We assume that the 
density of impurities is low enough that we can neglect 
impurity-impurity interactions and that there is no clus- 
tering. This is supported by experiments^^ on doped 
Bi2Te3 . Therefore, we add an on-site electron-impurity 
exchange interaction terni^ 



H = H, 



■ Si 



(4) 



where is the Hamiltonian for the electrons in a TI to 
be discussed in Sec. IV below. J is an exchange coupling 
constant, the sum extends over all impurity sites /, Sj 
and Sj are the spin operators of the impurity and electron 
on site /, respectively. 

We define the average magnetization of the impurities 
and electrons as M and m. The impurity magnetization 
is averaged over all impurities and over all sites. This vir- 
tual crystal approximation can be pictured as "spreading 
out" the localized magnetic moment of the impurities, so 
that there is a moment on all sites. We assume that 
the fiuctuations around the mean are small, and can be 
neglected, so that the interaction term is decoupled 



jj:sr- 



Si ~ JM 



Jm 



NJM ■ m(5) 



where A'^ is the number of sites. We assume that the 
magnetization is in the z direction, so M = Mz and 
m = mz. Then the mean-field Hamiltonian becomes, 
H^P = HfP + HfP - NJMm, where 



k 
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and 1-Le{k) is a matrix in spin and orbital spaces. Since 
the mean field Hamiltonian is decoupled, we can write 
the energies as a sum 



E{k, A/) = Ee{k) ^Ei- NJMm. 



(7) 



The first term refiects the energies of the mean field elec- 
tron Hamiltonian . The second term is a sum over 
the expectation values of the spin operator in the z di- 
rection at each impurity site Ei = Jm^jXi^ where A/ 
is the component of the spin in the z direction on im- 
purity /. The last term gives the overall shift in energy 
following from the mean-field decoupling Eq. ([5|. 

To find the magnetizations, we calculate the expec- 
tation values of the electron and impurity spins in the 
ensemble defined by ^ and find the equations 



M = -2SxBs{(3JmS), 

k,i 

where Bs{y) is the Brillouin functiorP^ 



(8) 



Bs{y) = coth y 



2S 



2S 



(9) 



U is the matrix that diagonalizes H^^ (it is a function of 
J, M, and k) and f{E) is the Fermi function. The chem- 
ical potential ji can be determined by summing over the 
average occupation for each energy (given by the Fermi- 
Dirac distribution) and equating to the total number of 
states. The number of states per site is n = f{Ee)/N. 
For example, for half filling (/i = 0) at T = for a model 
with four states per site, we have n = 2 since only half 
of those are occupied. The equations for m, M and fi 
are coupled non-linear equations which can be solved self 
consistently by an iterative procedure. 



B. Adding the surfaces 

Consider now a system with open boundary conditions 
in the z direction and periodic boundary conditions in 
the X and y directions (see Fig. [2|. The advantage of 
this setup is that in addition to the bulk, we can also 
investigate the magnetization near the two surfaces of 
the sample. We now take the real space Hamiltonian, 
and Fourier transform it only in x and keeping the 
dependence on z in real space. The problem is now a ID 
problem in as opposed to a OD problem before. 

Once again, we rewrite the interaction term and then 
define the magnetizations as 2;-dependent fields. The 
mean-field Hamiltonian is now = 
N^J^. Mj • mj, where j labels the layers and 



= E< [ne{k^)^JM 

Hf^ = J E • Si. 
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FIG. 2: A 3D TI in a slab geometry with a low density of 
randomly positioned magnetic impurities. 



Here 1-Le{k±) is a matrix in indices labeling the z coor- 
dinate in addition to spin and orbital spaces. M is a 
vector of length Lz (the number of sites in the z direc- 
tion), N± = N/Lz is the number of sites in the x-y plane, 
and hy M ^ (Jz we mean that M is expanded to a di- 
agonal matrix of dimension Lz and then we perform an 
external product with a^. We now diagonalize the full 
electron Hamiltonian He ) + /2 — /i, and find 
a matrix U which diagonalizes it. 

We calculate the electron and impurity magnetizations 
of each layer, which are given by the sum over the expec- 
tation values of the spin operators 



-2SxBs{(3JmiS), 



(11) 



where Uu is a d x d block of [/, and d is the dimension 
of the Hamiltonian in the bulk {d = 4 for Bi2Se3and 
d = 8 for perovskite, the models are defined below). The 
chemical potential can be determined as before, by solv- 
ing n = X]fc_L a f{^ociM'±))/N and a runs from 1 to Lzd. 



C. Estimate of the surface critical temperature 

Consider a ferromagnetically ordered 2D surface of a 
TI. We can make a rough estimate of the critical tem- 
perature. The effective Hamiltonian is given by ([T]) and 
([2|, and the energies are given by ([3|. The impurity 
and electron magnetizations are described by the cou- 
pled equations 



M = -2SxBs{f3JmS) 

f{Ek - /i) 



(12) 



NJ ^ dM 



For T j-surf haye m^M 0, so we expand the 
above equations to first order in M and m. We make use 
of dE/dM = J'^M/4:E and the assymptotic form of the 



FIG. 3: (a) A unit cell of the discretized Bi2Se3 model, a 
simple cubic lattice with two orbitals per site, (b) A unit 
cell of the perovskite lattice, which can be described as edge 
centered cubic. The four basis sites are labelled. 



Brillouin function Bs{y) — y{S + l)/35' for <C 1, and 
find 
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3^ 
JM 



M 



m 



-^S{S^l)/3Jxm, (13) 



We can evaluate the sum by converting it into an integral 
and imposing a momentum cutoff A. Assuming Ej^ c± 
±v\k\ (since M ^ at the critical point) we find a simple 
result-^^-^ for exact half filling, i.e /i = 0, 



nsurf 



3kB 



Aa\ J^x 
2^) XKv)' 



(14) 



For S 5/2, J 0.5eV, x 0.05, v = 2\^a and a 
cutoff Aa/(27r) 1/5 we find T^^^^^ lOOK. Away from 
half filling, when j3ii <^ 1^ there is a small correction 



-fi^/{4:kBvA) to the result in Eq. (14). 



A simple result can also be found for the case pji^ 1 
which is relevant for most values of ji inside the bulk gap. 



nsurf 



Ao 
27r 



J2 



(Aw)' 



:{^v-\^^\). (15) 



As expected, T^^^^ is seen to decrease away from half 
filling. 

To complete the argument we would now like to give 
a similar simple estimate for T^^^^. Unfortunately, the 
bulk critical temperature is not easy to estimate, since 
unlike the topologically protected surface state, whose 
physics is simple and universal, the bulk of a TI can be 
complicated and T^^^^ will generally depend on the de- 
tails of the band structure and other factors. For Bi2Se3 
with 5% concentration of Cr dopants, Ref. |32| estimates 
j-buik ^ YQj^ using first-principles numerical calculations. 
Comparing with our rough estimate for T^^^^ given above 
we see a clear indication that a T^^^^^ < T^""'^ regime can 
be easily obtained. 



IV. MODELS AND RESULTS 

Below we present our results for the two tight bind- 
ing models we considered in this study: a model for 
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FIG. 4: (a) Impurity spin expectation value as a function of temperature for the discretized Bi2Se3 lattice. We plot the order 
parameter on the surfaces ("surf"), the averaged bulk result for the 14 middle sites ("avg"), and the results of the separate 
bulk calculation ("bulk"), (b) Impurity spin expectation value as a function of the z coordinate, for different temperatures. 
For both plots the parameters are Lx = Ly = 40, Lz =20, 5* = 5/2, J = 0.5eV, x = 0.05, a nd the rest of the parameters 
are chosen to fit the Bi2Se3 dispersion close to the F point based on first principles calculations P^'^^^ 7o = 0.3391 eV,7± = 
0.0506eV,7^ = 0.0717eV,e = 1.6912eV,t± = 0.3892eV,t^ = 0.2072eV, A± = 0.2170eV and = 0.1240eV. 



Bi2Se3 regularized on a simple cubic latticd^and a model 
on the perovskite lattice^. These simple tight-binding 
models with spin-orbit coupling exhibit non-trivial topo- 
logical invariants for a broad range of model parameters 
and have been used widely in the literature to study the 
physical properties of topological insulators. 



Effective model for Bi2Se3 regularized on the 
cubic lattice 



We consider electrons hopping on a simple cubic lat- 
tice with two orbitals per site (Fig. The form 
of the Hamiltonian and the parameters (see caption 
Fig|4| were chosen to fit the dispersion near the cen- 
ter of the Brilloui n zone obtained by the first principles 
calculatioiPSmEl for Bi2Se3 . The electron Hamiltonian 
is given by 



ne{k) 



(16) 



do = rrik = e — 2^^^ti cos{kiai) 

i 

di = —2\i sin(/c^a^) 

= 70 - 2 ^ 7i cos{kiai) 

i 

where we have chosen the gamma matrices so that Fq = 
Ti (g) ctq, Fi = -Ts (g) 0-2, F2 = Ts (Ji, F3 = r2 (g ctq, and 
Ti are Pauli matrices in the orbital subspace and are 
Pauli matrices in the spin subspace. The energies are 



Ee = d^^^\ dl^dl 



dl^dl^ 5M 



(17) 



where M = JM/2, and 7, = ±1. Note that if there 
are no impurities ( J = 0) then we get the "clean" doubly 



degenerate electron spectrum Eq = d^^ y ^m* 

The results for the magnetization in this model are pre- 
sented in Fig[4j The bulk is ferromagnetically ordered up 
to T^^^^ :^ 73K, and the surface remains FM ordered up 
to T^""'^ 102K for the two surfaces (see Figji^). There- 
fore, the window in which the surface is ordered and the 
bulk is unordered (paramagnetic) is 29K. We plot the 
magnetization in the z direction in Fig[4]3. Here the effect 
can be seen clearly. If we ramp up the temperature from 
T = at first all spins, regardless of if they are in the 
bulk or surface, are fully polarized. As we increase the 
temperature, the magnetization of the bulk drops faster 
than the magnetization of the surface. Eventually we 
cross T^^^^ at which stage the magnetization of the bulk 
is zero, but the magnetization of the surface is finite. If 
we increase the temperature further, the magnetization 
of the surface drops, but the magnetization of the bulk 
remains at zero. Once we cross T^^^^, the magnetiza- 
tion of both the bulk and the surface vanish - thermal 
fluctuations have broken the ordered phases. 

In samples with open surfaces we observed spatial fluc- 
tuations in the bulk magnetization as the temperature 
approached T^^^^ from above, e.g. the T = 80K curve 
in Fig[4]3. We attribute these to the large magnetic sus- 
ceptibility of the bulk, which diverges at T ^ T^^^^, and 
the proximity of the ordered surfaces. The correct T^^^^ 
can be obtained from the bulk calculation with periodic 
boundary conditions, the results of which are plotted in 
Fig. [4^. 

We plot the critical temperature for the surface and the 
bulk as a function of the chemical potential in Fig [5] The 
surface critical temperature is maximal when the chem- 
ical potential intersects the surface Dirac point, which 
happens for /i 0.09eV, and falls on both sides, as ex- 
pected. The result agrees well with the predicted linear 
dependence of the surface critical temperature on ja pre- 
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FIG. 5: Critical temperature as a function of the chemical 
potential /i for the discretized Bi2Se3 lattice. We plot the sur- 
face critical temperature (squares, blue), and the approximate 
surface temperature from Eq. (15) (solid line, black). In ad- 



dition, we plot the bulk critical temperature as obtained from 
an average over the bulk sites (circles, red). The error bars 
were obtained by fitting the magnetization below the critical 
temperature to the known behaviour for a second order phase 
transition, M oc (Tc — T)^/^. The parameters are as in Fig. [4] 
except for = 10. The cutoff was chosen for best fitting by 
eye Aa/27r = 0.215. The vertical dashed lines mark the bulk 



dieted in Eq. (15). The bulk critical temperature is ap- 



proximately constant within the bulk gap, as expected 
- since there are no bulk states within the gap, chang- 
ing the chemical potential should not change the critical 
temperature. 



B. Perovskite lattice model 

To ascertain whether the results obtained for the 
Bi2 Sea model are in fact generic, we investigate another 
simple lattice model known to give a robust TL The per- 
ovskite lattice, or edge-centered cubic lattice, consists of 
atoms on a simple cubic lattice with a four point ba- 
sis, corresponding to atoms at the center of each edge 
(Fig. [sJd). Tight binding electrons on this lattice (with 
spin orbit coupling) were recently shown to form a TI.^^ 
The momentum-space electron Hamiltonian, written in 
the basis of four sublattice sites indicated in Fig. |3]3, is 
T-Le = ^0 ^SO' Here Ho is the hopping part 



no = -2t 



I cos/ca; cos/cjy cos/c^ 
COS kx ^ 
cos/c^ 

VcosA:^ 



(18) 



and 1-iso is the spin orbit part, 1-Lso = ^ 
and 



Hso — 4Ai sin ky sin kz 



/O ON 



1 

Vo -1 Oy 



50' 



(19) 



and similarly for y and z. 

The results for magnetization in this model are pre- 
sented in Fig|6] In this case we observe that for a large 
range of temperatures, the magnetization on basis sites 
1 and 4 is similar and opposite in sign from the mag- 
netization on basis sites 2 and 3 (which are equal due 
to a rotational symmetry around the z axis). This mo- 
tivates the definition of a FM order parameter Sj^ = 
(6'i^i+5'i^2 + 'S'i^3 + 6'i^4)/4 and an antiferromagnetic (AF) 
order parameter Sf^^ = (5'i,i — 5'i,2 — 'S'i,3 + 6'i,4)/4 (where 
Si^i is the expectation value of the impurity spin on ba- 
sis site ^ of lattice site i). We see that the system is 
ferromagnetically ordered up to T 4.5K, where it be- 
comes AF ordered. The bulk remains AF ordered until 
j-buik ^ 50. 5K, and the surface remains AF ordered until 
j^surf _ and 78K for the two surfaces. Therefore, the 
maximum window in which the surface is ordered and the 
bulk is unordered (paramagnetic) is c:^ 17. 5K. Note that 
the difference in surfaces results from the fact that the 
unit cell is not symmetric under reflection along z - the 
"top" of the system ends with basis site 4 and the bot- 
tom ends with basis site 1 (see Fig. 3|i). In addition, we 
see that the surfaces undergo an additional partial phase 
change signified by discrete jumps of the order parame- 
ter seen at T 4.5, 5.5, 13.5, 35. 5K, which we claim as 
additional evidence that the bulk and surface differ. 



V. CONCLUSIONS 

We have demonstrated that magnetically doped topo- 
logical insulators can have a sizeable window where the 
bulk is paramagnetic and the surface is magnetically or- 
dered. Our conclusions are based on general arguments 
that involve only universal properties of the topologically 
protected surface states and on numerical calculations 
performed on simple lattice models of 3D topological in- 
sulators. Physically, these results are in accord with the 
intuition that the metallic state on the surface should 
be more susceptible to magnetic ordering than the insu- 
lating bulk, although this expectation is only partially 
borne out in systems with strong spin-orbit coupling. 

The results reported in this study rely on two key 
approximations: (i) the mean- field decoupling of the 
exchange interaction between electrons and magnetic 
dopants indicated in Eq. ([5|, and (ii) the 'virtual crystal' 
approximation which replaces the localized magnetic mo- 
ments of dopant atoms by their average over all lattice 
sites. We have attempted to bypass the latter approxi- 
mation by solving the problem in real space without aver- 
aging over sites. This procedure gives reasonable results 
for T^buik compared to those obtained in the virtual crys- 
tal but is numerically more costly. The system sizes that 
we could simulate did not allow us to unambiguously de- 
termine the surface critical temperature. Nevertheless 
based on these results we feel that, in a large crystal, (ii) 
will not lead to a significant error in the determination 
of the critical temperatures. 



7 




FIG. 6: (a) Ferromagnetic order parameter as a function of temperature for the perovskite lattice, with parameters Lx — Ly — 
20, Lz = 10, 5* = 5/2, J = 0.25eV, t = leV, A = 0.2eV, x = 0.05. We plot the order parameter on each of the surfaces 
("surfl" and "surf2"), the result for the site dX z — Lz/2 ("z"), and the results of the separate bulk calculation ("bulk"), (b) 
Antiferromagnetic order parameter as a function of temperature for the perovskite lattice, with parameters as in (a). Here we 
also plot the average over the four central sites ("avg"). 



Since fluctuations around the mean field result are 
typically stronger in 2D than in 3D, they will likely re- 
duce the size of the temperature window between T^^^^ 
and T^""'^. On the other hand, our mean field calcula- 
tion did not include electron-electron interactions, which 
would tend to stabilize the long range magnetic order and 
thereby strengthen the ordered phases. It has been shown 
recently that electron-electron interactions can lead to 
spontaneous breaking of TRI on the surface of a TI,'^ 
even in the absence of magnetic dopants. Thus the inter- 
actions might in fact strengthen the effect found in our 
study, although a more detailed investigation would be 
needed in order to obtain a quantitative result. Overall, 
in our opinion it is very likely that the combined effect of 
the magnetic dopants and the electron-electron interac- 
tions can account for the experimentally observed surface 
excitation gap without bulk magnetic order in Mn and 
Fe doped Bi2Se3 . 

We note that the critical temperature for the bulk 
magnetic ordering of the magnetically doped topologi- 
cal insulator Bi2-ajMna;Te3 was recently measured^^ to 
be 9 - 12K for x = 0.04 - 0.09. This critical temper- 
ature is smaller than our results and those of Ref. \32\ 
which estimates T^^^^ c:^ 70K for Cr-doped Bi2Se3 using 



first-principles numerical calculations. We chose our cou- 
pling constant J based on the first-principles results, so 
we expect that for a reduced coupling constant the tem- 
perature window would shift to lower temperatures. As 
noted above, our arguments are universal, hence the ex- 
act details of the material and coupling might alter the 
result quantitatively but should not change it qualita- 
tively. Taking the above finding for Bi2-a;Mna;Te3 as a 
guideline, one may surmise that T^^^^ for Fe and Mn 
doped Bi2Se3 lies in a similar range of temperatures. It is 
then entirely possible that the ARPES experiments j ^^ * ^^ ^ 
performed at ~ 20K, detected surface magnetic order- 
ing without bulk magnetism as advocated in this paper. 
Careful surface-sensitive magnetic measurements, as pro- 
posed recently,^^ might be able to probe this intriguing 
phenomenon directly. 
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